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FIGURES 


(a)  Typical  alkali  ground  state  hyperflne  transition  llneshape 
in  the  presence  of  a  few  Torr  of  buffer  gas.  (b)  Full  Doppler 
profile  of  the  alkali  ground  state  hyperflne  transition,  observed 
in  the  absence  of  buffer  gas.... . . . 

Linewidth  of  ®^Rb  ground  state  hyperflne  transition  as  a 
function  of  the  ratio  of  a,  the  atom's  mean  free  path,  to  X, 
the  transition  wavelength . 


Hyperflne  transition  llneshape  resulting  from  trajectory 
analysis  of  a  Kb  atom  in  a  wall-coated,  bufferless  cell 
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(a)  Hyperflne  llneshape  for  a  °  Rb  atom  in  a  0.34-cm-radius, 
wall-coated,  bufferless  cell  obtained  from  the  trajectory 
calculation,  (b)  ®'Rb  hyperflne  llneshape  after  relationship 
between  precollision  velocities  and  postcolllslon  velocities, 
found  in  the  wall-coated,  bufferless  cell,  has  been  ignored... 


Hyperflne  lineshapes  for  ®^Rb  atom  in  cells  of  varying  radii. 


1 .  INTRODUCTION 


The  phenomenon  of  producing  sub-Doppler  spectral  linewidths  in  atomic 
emissions  by  confining  atoms  with  a  buffer  gas  is  well  known.  This  line¬ 
narrowing  effect  was  first  analyzed  by  R.  H.  Dicke^  and  is  typically  termed 
Dicke  narrowing.  Dicke  considered  atoms  confined  by  a  nonperturbing  buffer  gas 

of  sufficiently  high  pressure  that  the  mean  free  path  of  the  emitting  atom  was 

2_ 

much  less  than  the  wavelength  of  the  emitted  radiation.  Subsequent  analyses 

have  relieved  this  constraint  and  have  produced  a  formalism  that  predicts  a 

Doppler-broadened  Gaussian  lineshape  in  the  limit  of  low  buffer  gas  pressure 

and  that  reproduces  Dicke 's  result:  a  narrowed  Lorentzian  lineshape  in  the 

limit  of  high  buffer  gas  pressure.  One  of  the  most  dramatic  realizations  of 

Dicke  narrowing  is  observed  in  the  magnetic  dipole  transitions  of  the  hyper- 

87 

fine  levels  of  alkali  atoms.  The  ground  state  hyperfine  transitions  of  Rb 

1  OO 

and  Cs,  at  frequencies  of  6.8  GHz  and  9.2  GHz,  respectively,  have  been 
extensively  studied  because  of  their  applications  in  atomic  frequency 
standards.7-^  In  Fig.  1  we  compare  a  typical  lineshape  that  might  be 
observed  for  the  ground  state  hyperfine  transition  in  °  Rb  in  the  presence  of 
a  few  Torr  of  buffer,  with  a  calculation  of  the  full  Doppler  profile. 

Typically  observed  linewidths  are  about  300  Hz,  and  only  a  small  fraction  of 
that  due  to  the  residual  Doppler  contribution;  the  full  Doppler  width  is 
computed  to  be  around  10  kHz. 

The  lineshape  derived  by  Dicke^  is  given  by 

Km)  _  2 it  D/X2 _ 

Xo  (w-w,)2  +  (2#  D/X2)2 

where  D  is  the  diffusion  constant  of  the  atom  through  the  particular  buffer 

gas,  and  X  is  the  wavelength  of  the  transition,  X  *  c/2ttuj  .  The  full  width  at 

2  Q7  ° 

half  intensity  of  the  line  is  4nD/X  (rad/s).  For  Rb  in  a  typical  noble 
gas  buffer,  e.g.,  Kr,  one  computes  the  residual  Doppler  width,  i.e.,  the 
Dicke-narrowed  contribution,  to  be  approximately  40  Hz/p,  where  p  is  the 


pressure  of  the  buffer  In  Torr.  In  a  few  Torr  of  buffer  gas,  then,  the  Dicke- 
narrowed  contribution  is  only  a  very  small  fraction  of  the  total  linewidth. 

In  the  usual  experiments,  the  major  contributions  to  the  linewidth  arise  from 
the  effects  of  the  optical  pumping  radiation,  the  microwave  field,  and 
magnetic  field  gradients  and  the  like. 

It  is  also  well  known  that  sub-Doppler  linewidths  can  be  obtained  in 
absorption  cells  containing  no  buffer  gases  but  having  wall  coatings,  such  as 
paraffin,  to  minimize  the  relaxation  effects  of  wall  collisions  on  the  atomic 
spin  orientation. Microwave  linewidths  as  narrow  as  11  Hz  have  been 
observed  in  some  of  these  experiments.1  The  spectral  narrowing  of  the  rf  or 
microwave  lines  in  these  buffer less,  wall-coated  cells  has  been  explained  by 
appealing  to  the  example  of  an  emitting  atom  confined  to  a  one-dimensional 
box.  In  this  model,  an  atom  is  confined  to  a  box  of  dimension  a,  travels  at 
constant  velocity  V,  and  is  allowed  only  to  rebound  elastically,  reversing 
direction  upon  a  wall  collision.  It  can  be  easily  shown  that  for  nonintegral 
values  of  a/A  the  atom  emits  at  both  the  normal  unshifted  frequency  and  the 
sideband  frequencies  determined  by  its  velocity.  For  an  ensemble  of  atoms 
with  a  Maxwell-Boltzmann  velocity  distribution  confined  to  a  one-dimensional 
box,  one  finds  the  llneshape  to  be  composed  of  a  central  spike  at  the  normal 
unshlfted  frequency,  and  a  broad  pedestal  built  up  from  the  velocity-dependent 
sideband  frequencies. 

In  reality,  however,  atoms  confined  to  bufferless,  wall-coated  cells  are 
not  well  described  by  the  one-dimensional  box  model.  All  real  cells  are  three 
dimensional.  Simply  generalizing  the  one-dimensional  model  to  two  dimensions 
changes  the  physical  problem  significantly.  First,  the  distance  that  the  atom 
may  travel  between  collisions  is  no  longer  fixed;  rather,  it  is  distributed 
from  zero  to  the  diameter  of  the  cell.  Then,  even  if  the  atom  is  assumed  to 
be  emitted  from  the  cell  surface  with  a  fixed  speed  S,  a  range  of  velocity 
components  is  still  possible  as  the  atom  rebounds  randomly  within  the  cell. 

In  two  dimensions  the  well-defined  sideband  emission  frequencies  disappear. 

The  range  of  velocities  causes  sideband  frequency  power  to  be  spread  over  all 
frequencies  between  *  Su>Q/2nc.  Averaging  these  lineshapes  over  a  Maxwell- 
Boltzmann  distribution  results  in  a  center  peak  resting  on  a  non-Doppler 


pedestal.  Furthermore,  purely  elastic  collisions  with  cell  walls  that  allow 
atoms  to  bounce  with  no  change  in  speed  are  not  physically  realistic.  In 
fact,  coated  cell  walls  tend  to  be  somewhat  sticky.  Atoms  that  collide  with 
the  coated  wall  stick  for  some  period  of  time,  between  10“**  and  10“*®  s>13-14 
and  then  are  re-emitted  at  some  random  angle  and  random  velocity  governed  by 
the  respective  angle  and  velocity  distribution  functions.  Relaxing  the 
condition  of  elastic  collisions  at  the  wall  in  the  one-dimensional  model  also 
results  in  a  sharp  central  spike  resting  on  a  clearly  non-Doppler  pedestal 
because,  once  again,  each  atom  experiences  a  number  of  different  velocities 
over  its  lifetime  and  thus  would  emit  at  numerous  sideband  frequencies  so  that 
the  spectral  power  is  spread  out,  leaving  the  greatest  intensity  at  the  normal 
unshifted  frequency.  For  the  various  reasons  outlined  the  "particle  in  a  box" 
analysis  of  motional  narrowing  in  bufferless,  wall-coated  cells  is 
unsatisfying. 

One  aspect  of  Dicke's  criterion  for  collisional  lineshape  narrowing*  is 
that  the  emitting  atom  suffers  collisions  that  do  not  perturb  the  atom's 
internal  state.  As  indicated  by  Vanier  et  al.,*®  collisions  of  emitting  atoms 
with  a  buffer  gas  may  be  described  in  the  same  manner  as  collisions  with  a 
wall  surface,  both  satisfying  Dicke's  criterion.  Thus  an  atom  undergoing 
hyperfine  transitions  should  not  really  "know"  if  it  is  confined  by  a  buffer 
gas  or  a  wall,  and  consequently,  one  should  not  need  two  apparently  distinct 
formalisms,  buffer  gas  and  particle  in  a  box,  to  describe  spectral  line 
narrowing  under  collisional  confinement.  However,  a  naive  application  of  the 
line-narrowing  formalism  to  the  bufferless,  wall-coated  cell  completely 
fails.  In  Fig.  2  we  plot  the  full  width  at  half  intensity  of  the  residual 
Doppler  contribution  to  the  °  Rb  microwave  lineshape  at  6.8  GHz,  as  predicted 
by  the  complete  line-narrowing  formalisms,  versus  the  dimensionless  parameter 
a/A,  where  a  is  the  mean  free  path  of  the  atoms.  The  area  to  the  left  of  the 
dashed  line  in  the  figure  indicates  the  region  of  applicability  of  Dicke's 
result.  To  the  right  of  the  dashed  line,  the  linewidth  is  seen  to  smoothly  go 
to  the  full  Doppler  width,  as  expected.  In  addition,  we  have  also  plotted  in 
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Figure  2.  Linewidth  of 


Linewidth  of  °  Rb  ground  state  hyperfine  transition  as  a  function 
of  the  ratio  of  a,  the  atom's  mean  free  path,  to  X,  the  transition 
wavelength.  Linewidths  were  calculated  according  to  Galatry.  Th 
dashed  line  represents  the  limit  of  validity  of  Dicke's  original 
analysis.*  The  "+"  symbol  corresponds  to  the  estimate  of  the 
residual  Doppler  contribution  to  the  measured  linewidth  in  Ref.  12 


Fig.  2  the  estimate  of  the  residual  Doppler  contribution  to  the  linewidth 
observed  by  Robinson  and  Johnson,  denoted  by  the  symbol  In  this 

particular  case  we  took  the  mean  free  path  of  the  Rb  atoms  to  be  the  radius  of 
the  absorption  cell.  It  is  quite  apparent  that  the  standard  formulations  of 
line  narrowing  fail  dramatically  in  this  situation. 

We  have  investigated  the  line-narrowing  phenomenon  in  a  bufferless,  wall- 
coated  absorption  cell  following  the  same  formalisms  developed  for  the  buf¬ 
fered  cell.  Particle  trajectories  within  the  absorption  cell  were  generated 
using  Monte  Carlo  techniques.  This  method  of  analyzing  atomic  motion  in  a 
spherical  cell  has  been  applied  previously  to  motional  narrowing  of  Zeeman 
lineshapes  in  the  presence  of  an  inhomogeneous  magnetic  field.  *  This  form 
of  analysis  is  particularly  appropriate  for  hyperfine  transitions,  because 
wall  coatings  exist  that  minimize  hyperfine  state  changes  during  the  wall 
collision.  We  have  found  that,  although  the  trajectories  followed  by  the  atom 
over  its  hyperfine  state  lifetime  have  a  statistical  nature  because  atoms 
striking  the  coated  wall  are  adsorbed  for  a  finite  time  and  then  re-emltted 
randomly,  there  still  is  a  very  high  probability  that  the  Doppler  shifts 
before  and  after  a  particular  wall  collision  are  of  opposite  sign.  This  leads 
to  an  effective  cancellation  of  Doppler-shif t  effects.  The  resulting 
lineshape  displays  a  sharp  central  spike  resting  on  a  clearly  non-Gaussian 
pedestal.  Additionally,  analysis  of  the  atomic  trajectories  in  a  wall-coated 
cell  demonstrates  a  correlation  between  postcollision  velocities  and  inter¬ 
collision  times.  While  not  of  great  significance  in  these  cells,  a  similar 
|  correlation,  albeit  a  much  weaker  one,  exists  in  the  buffered  cell.  This 

correlation,  which  has  been  completely  neglected  in  previous  analyses  of  the 
line-narrowing  phenomenon,  will  lead  to  additional  narrowing  of  the  residual 
Doppler  contribution  beyond  that  given  in  the  standard  analyses.  A  major 
|  value  of  our  new  analysis  Is  that  spectral  line  narrowing  due  to  confining 

collisions  can  be  explained  by  a  single  formalism  no  matter  what  the  origin  of 
the  confining  collisions. 
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2.  ANALYSIS  AND  RESULTS 


An  effective  method  of  analyzing  certain  aspects  of  radiation  problems, 
as  pointed  out  by  Weisskopf,  is  to  treat  each  radiating  atom  as  a  classical 
oscillator.  This  technique  has  been  exploited  in  previous  analyses  of  line 
narrowing  for  an  atom  surrounded  by  a  buffer  gas.  The  displacement  of  the 
oscillator,  x(t),  is  given  by 

oi  t 

x(t)  =  A  cos[u>  t  +  —  /  V  (t)  dt]  (2.1) 

OCX 

o 

with  <»>  the  atomic  transition  frequency,  Vv  the  velocity  component  leading  to 
the  Doppler  shift,  c  the  speed  of  light,  and  A  a  constant.  For  the  time  T(i) 
between  the  (i-l)th  and  ith  collisions  with  the  buffer  gas  atoms,  the  radiat¬ 
ing  atom  has  a  fixed  velocity  Vx(i).  This  allows  Eq.  (2.1)  to  be  rewritten  as 

oj  N 

x(t)  *  A  cos  [w^t  +  —  E  Vx(i)  T(i)]  (2.2) 

with  N  being  the  number  of  collisions  in  time  t. 

Typically,  the  postcollision  velocities  and  intercollision  times  are 
treated  as  independent,  stationary,  random  variables.  The  autocorrelation 
function  of  the  oscillator,  F(x),  is  given  by^ 

F(t)  =  <  x(o)*x(T)>  (2.3) 

v. 

where  the  brackets  indicate  stochastic  averaging  over  the  random  variables  Vx 
and  T.  To  perform  the  averaging  requires  the  probability  density  functions 
for  the  random  variables.  For  the  radiating  atom  surrounded  by  a  buffer  gas 
the  postcollision  velocities  are  distributed  according  to  the  Maxwell- 
Boltzmann  distribution;  the  intercollision  times  are  assumed  to  be  distributed 
according  to  a  simple  Poisson  distribution.^  In  this  model,  the  velocity,  and 
hence  the  Doppler  shift  after  a  collision,  is  assumed  to  be  completely 
independent  of  those  prior  to  the  collision.  The  Fourier  transform  of  the 


autocorrelation  function  Is  the  oscillator's  power  spectral  density  or,  in 
this  case,  would  be  more  appropriately  termed  the  transition  lineshape.  Line- 
widths  obtained  using  this  procedure  yield  the  solid  curve  plotted  on  Fig.  2. 

2.1  ANALYSIS  OF  MOTIONAL  NARROWING  IN  THE  WALL-COATED  CELL 

It  has  been  shown  that  the  naive  application  of  the  results  of  the  stand¬ 
ard  buffer  gas  analysis  to  motional  narrowing  in  the  wall-coated  cell  does  not 
yield  the  correct  linewidth.  Carrying  out  a  more  detailed  analysis  of  the 
wall-coated  cell  case  within  the  stochastic  averaging  framework  is  not 
straightforward,  because  (1)  while  the  postcollision  velocities  will  undoubt¬ 
edly  be  distributed  in  a  Maxwell-Boltzmann  form,  it  is  not  obvious  that  the 
intercollision  times  are  still  distributed  according  to  a  Poisson  distribu¬ 
tion,  and  (2)  the  assumption  that  there  is  no  relation  between  the  velocity 
and  Doppler  shift  before  and  after  a  collision  is  no  longer  even  approximately 
correct.  In  a  coated  cell  there  is  a  high  probability  that  the  velocities 
before  and  after  a  collision  will  have  opposite  signs  and  their  Doppler  shifts 
will  tend  to  cancel.  In  effect,  a  Markoff ian  relation  exists  between  the 
velocities  before  and  after  a  collision.  Aspects  of  the  effects  of  similar 
relationships  have  been  discussed  by  Anderson^  and  Kubo^  with  regard  to 
motional  narrowing  in  nuclear  magnetic  resonance  (NMR)  spectra.  However, 
attempting  to  apply  their  analyses  to  this  problem  is  quite  difficult. 

An  alternative  to  these  standard  stochastic  analyses  is  to  explicitly 
obtain  x(t)  as  a  function  of  time  and  then  perform  its  Fourier  analysis.  To 
carry  this  out,  a  series  of  Vx(i)'s  and  T(i)'s  as  would  be  observed  in  a  wall- 
coated  cell  must  be  generated.  This  has  been  done  by  calculating  the 
classical  trajectories  of  a  particle  in  the  wall-coated  cell.  A  similar 
method  of  analysis  was  applied  by  Watanabe  et  al.  in  their  study  of  the 
motional  narrowing  of  Zeeman  lineshapes  for  atoms  in  inhomogeneous  magnetic 
fields.  They  found  this  Monte  Carlo  analysis  to  be  quite  effective  for  their 
problem. 

It  is  well  known  that,  in  collisions  between  alkali  atoms  and  cell  walls 
coated  with  materials  like  paraffin,  the  colliding  atoms  are  physically 
adsorbed  and  then  re-emitted.^’ ^  The  direction  of  re-emission  is 


characterized  by  two  angles,  9  and  $.  Angle  0  ranges  from  -n/2  to  w/2  with 
respect  to  the  direction  normal  to  the  cell  wall  at  the  point  the  particle  has 
impacted.  Angle  $  is  the  azimuthal  angle  about  the  normal  direction  and 
varies  from  0  to  2n.  The  angle  9  has  been  found  to  be  cosine  distributed;  *  J 
♦  is  uniformly  distributed  throughout  its  range.  The  cos  (9)  dependence 
results  in  equal  probabilities  for  traveling  from  the  initial  surface  point  to 
any  other  point  on  the  sphere  surface.  Consequently,  given  the  initial  point, 
the  next  impact  may  be  determined  by  randomly  selecting  a  point  on  the  surface 
of  the  sphere.  The  emission  speed  S  was  assumed  to  satisfy  the  Maxweli- 
Boltzmann  speed  distribution.  Generation  of  a  random  variable  satisfying  this 
distribution  was  performed  using  standard  techniques.  »  Once  a  particular 
S  and  the  next  impact  point  are  selected,  it  is  possible  to  determine  the 
intercollision  time  and  the  component  of  velocity  in  the  direction  of  the 
detector,  Vx(i),  that  leads  to  the  Doppler  shift.  Repeating  this  procedure 
enables  generation  of  a  full  particle  trajectory. 

To  account  for  the  finite  coherence  lifetime  tQ  of  the  radiating  atom, 

Eq.  (2.2)  assumes  the  form 
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*Recent  experimental  measurements  performed  by  L.  R.  Martin  and  C.  M.  Kahla  of 
The  Aerospace  Corporation  show  that  Cs  particles  recoiling  from  a  paraffin 
surface  satisfy  a  cos  (9)  distribution.  A  similar  emission  distribution 
should  be  satisfied  by  Rb  atoms  recoiling  from  a  paraffin  coating  in  the 
standard  wall-coated  cell. 


where  Che  Vx(i)'s  and  T(i)'s  are  from  Che  CrajecCory  calculation,  Eq.  (2.2). 
Racher  chan  calculaCing  Che  auCocorrelacion  funecion  of  x(c)  and  Chen 
obcaining  Che  lineshape  from  ics  Fourier  transform,  we  performed  Che  direcC 
Fourier  analysis  of  x(c).  The  lineshape  G(u> )  is  relaced  direcdy  Co  x(c)  in 
Eq.  (2.5) :28 

G(a>)  -  |  /  e_ia)t  x(c)  dc|2  (2.5) 

The  Fourier  Cransform  of  x(c)  was  compuCed  numerically  and  was  facilicaced  by 
nocing  chac  beCween  collisions  Che  radiacion  from  Che  aeon  appears  Co  have  a 
fixed  frequency  and  ChaC  Che  frequency  changes  discreCely  upon  each  colli¬ 
sion.  Equacion  (2.5)  can  Chen  be  broken  inCo  a  sum  of  inCegrals,  one  for  each 
inCercollision  period,  wich  each  of  Che  inCegrals  capable  of  being  evaluaCed 
analycically. 

SCricCly  speaking  one  should  expliciCly  Cake  inCo  accounC  Che  face  chac 
CQ  represenCs  an  average  llfecitne  of  Che  aCom  and  weighc  Che  collisional 
effects  over  Che  acorn's  lifecime  using  an  exponential  funecion.  Because  Che 
average  lifecime  is  so  much  larger  Chan  Che  average  inCercollision  cine  in  Che 
syscems  of  inCeresC,  only  a  very  few  aComs  do  noc  parcicipace  in  Che  colli¬ 
sional  averaging.  As  will  be  shown,  collisional  narrowing  in  chese  systems  is 
extremely  effective,  and  thus  (for  clarity)  we  decided  to  consider  only  an 
average  atomic  state  lifetime,  knowing  Chat  under  the  specified  conditions  the 
error  will  be  negligible.  However,  in  a  more  general  situation  in  which  a 
substantial  fraction  of  the  atoms  have  relatively  few  collisions  over  their 
lifetimes,  a  properly  weighted  lifecime  funecion  would  have  to  he  incorporated 
into  Eq.  (2.4)  and  appropriate  averaging  performed.  One  further  computational 
aspect  should  also  be  discussed.  Oscillator  trajectory  x(t)  is  specified  by  a 
trajectory  lasting  a  period  tQ.  In  principle,  an  infinite  number  of  different 
trajectories,  each  lasting  tQ,  can  be  generated.  Each  trajectory  will  yield  a 
different  x(t)  with  a  slightly  different  G(iv).  To  obtain  the  desired 
accuracy  in  the  calculated  lineshapes ,  several  hundred  trajectories  were 
averaged  together,  corresponding  to  approximately  100,000  wall  collisions  for 
each  lineshape. 
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Application  of  this  computational  procedure  to  °'Rb  at  25°C  in  a  wall- 
coated  spherical  cell  with  a  radius  of  3.4  cm  yields  the  llneshape  shown  in 
Fig.  3.  The  linewidth  Ls  approximately  10  Hz,  indicating  that  in  the  wall- 
coated  cell  Dicke  narrowing  effectively  removes  all  Doppler  broadening  and 
yields  the  linewidth  specified  by  the  radiating  atom's  coherence  lifetime  in 
the  cell,  approximately  0.1  s.  This  result  is  in  agreement  with  the  measure- 
ments  of  Robinson  and  Johnson,  in  which  no  residual  Doppler  broadening  was 
observed.  The  small  oscillations  in  the  figure  are  due  to  the  computational 
model,  in  which  the  atomic  radiation  is  stopped  abruptly  at  tQ.  Averaging 
llneshapes  of  this  type,  generated  from  an  exponential  distribution  of  radia¬ 
tive  lifetimes,  would  yield  the  expected  Lorentzian  llneshape.  Additionally 
the  central  peak  is  found  to  be  resting  on  a  broad,  non-Gaussian,  pedestal 
with  a  width  of  approximately  two-thirds  that  of  the  calculated  Doppler 
width.  Published  experimental  observations  of  pedestals  are  limited;  however, 
Robinson  and  Johnson10  report  a  Gaussian  pedestal,  whose  width  equalled  the 
full  Doppler  linewidth.  Although  the  trajectory  formalism  that  we  have 
described  provides  an  excellent  model  for  linewidth  narrowing,  there  is  at 
present  some  difficulty  in  reproducing  the  pedestal  portion  of  the  line  as 
experimentally  observed.  We  discuss  the  implications  of  this  problem  in 
Section  2.3. 

The  extremely  efficient  Dicke  narrowing  yielding  the  central  peak  found 
in  the  wall-coated,  bufferless  ceil  is  due  primarily  to  the  high  probability 
that  the  velocity  after  a  collision  will  have  the  opposite  sign  to  that  prior 
to  the  collision,  thereby  cancelling  the  Doppler  shifts.  The  significance  of 
this  relationship  between  the  pre-  and  postcoilision  velocities  may  be  readily 
demonstrated  by  repeating  the  trajectory  calculation.  However,  this  time, 
upon  each  collision,  the  position  of  the  particles  is  randomized  over  cell 
surface.  While  each  trajectory  is  still  subject  to  the  spatial  constraints  of 
the  ceil  and  the  cosine  emission  distribution,  any  "memory"  of  a  prior  vel¬ 
ocity  is  removed. 

Figure  4a  depicts  the  llneshape  for  a  particle  in  a  0.34-cm-radius ,  wall- 
coated,  bufferless  cell  obtained  from  the  standard  trajectory  calculation. 
Beneath  it,  Fig.  4b  indicates  the  llneshape  resulting  from  the  modified 
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Figure  A.  (a)  Hyperfine  llneshape  for  a  ®^Rb  atom  in  a  0.34-cra-radius ,  wall- 
coated,  bufferless  cell  obtained  from  the  trajectory  calculation. 
The  linewidth  is  10  Hz.  (b)  ®'Rb  hyperfine  lineshape  after  rela¬ 
tionship  between  precollision  velocities  and  postcollision  veloc¬ 
ities,  found  in  the  wall-coated,  bufferless  cell,  has  been  ignored 
The  linewidth  is  approximately  1400  Hz. 


trajectory  calculation  just  described.  The  0.34-cm  radius  was  chosen  for 
clarity  of  example  such  that  the  uncorrelated  linewidth  would  still  be  well 
below  the  full  Doppler  width.  The  linewidth  in  Fig.  4b  is  approximately 
1400  Hz,  over  two  orders  of  magnitude  larger  than  the  10-Hz  linewidth  of 
Fig.  4a.  The  linewidth  displayed  in  Fig.  4b  is  essentially  the  same  as 
predicted  by  the  standard  statistical  analyses  of  Dicke  narrowing  in  a  buffer 
gas.  Subtle  differences  between  the  lineshape  displayed  in  Fig.  4b  and  that 
given  by  the  standard  analyses  are  due  to  the  non-Poisson  intercollision  time 
distribution  in  the  wall-coated  cell  and  thus  an  exact  comparison  with  Fig.  2 
is  not  appropriate.  This  calculation  clearly  verifies  that  it  is  the  rela¬ 
tionship  between  the  precollision  and  postcollision  velocities  that  leads  to 
the  extremely  effective  narrowing  in  a  wall-coated,  bufferless  cell. 

2.2  ANALYSIS  OF  MOTIONAL  NARROWING  IN  THE  BUFFER-GAS  CASE 

Examination  of  the  trajectories  of  atoms  in  the  wall-coated,  bufferless 
cell  shows  that  there  is  a  correlation  between  the  postcollision  velocity  of 
an  atom  and  the  intercollision  time.  Faster-moving  atoms  typically  have 
shorter  intercollision  times.  For  the  bufferless  cell,  the  trajectory 
analysis  includes  these  effects.  Any  similar  correlation  in  the  buffer  gas 
cell  has  typically  not  been  taken  into  account  in  the  standard  Dicke 
calculations.  The  standard  statistical  analysis  of  motional  narrowing  for  a 
particle  surrounded  by  a  buffer  gas  assumes  that  the  postcollision  velocities 
and  intercollision  times  are  completely  independent.  While  this  is  a  fairly 
good  assumption  it  is  not  rigorously  correct.  Some  correlation  exists  even  i 
the  buffer-gas  case,  albeit  not  nearly  as  strong  as  that  in  the  wall-coated 
cell.  Clearly,  the  velocity  of  the  radiating  atom  will  affect  to  some  extent 
the  intercollision  time. 

Typically,  intercollision  times  are  assumed  to  obey  the  simple  Poisson 
distribution,  and  the  most  probable  Intercollision  time  is  the  inverse  of  the 
average  collision  frequency.  Most  probable  intercollision  time  Tp  does  not 
depend  on  the  precise  speed  of  the  radiating  atom  after  each  collision  but  is 
a  function  of  the  speed  of  the  radiating  atom.  The  relation  between  the  most 
probable  collision  time  of  a  radiating  atom  with  a  precise  speed  S  surrounded 


by  a  buffer  gas  of  mass  m,  at  temperature  T,  and  with  a  number  density  N  is 
given  by^® 

(ip)'1  -  N  tr a1  (if 8)_l/2  (81/2S)_1  i|»  (01/2S)  (2.6) 
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and  k  the  Boltzmann  constant. 

To  analyze  the  effects  of  correlation  in  the  buffer-gas  case,  x(t)  as 

specified  in  Eq.  (2.4)  is  generated.  After  each  collision  a  random  speed  S(i) 

is  generated  by  taking  the  magnitude  of  a  velocity  vector  whose  three 

components  were  randomly  selected  from  independent  Maxwell-Boltzmann  velocity 

distributions.  This  clearly  is  a  simplistic  model,  considering  that  the 

postcollision  speed  may  be  strongly  related  to  the  precollision  speed. 

However,  it  is  felt  to  be  adequate,  because  our  primary  interest  is  in  the 

relationship  between  postcollision  speed  and  intercollision  time — two 

quantities  typically  taken  to  be  uncorrelated.  Velocity  component  Vx(i) 

generates  the  Doppler  shift  after  the  ith  collision.  The  Intercollision  times 

are  chosen  as  follows.  Given  the  postcollision  speed,  a  most  probable 

intercollision  time  r  (i)  is  calculated  using  Eqs.  (2.6)  through  (2.8).  This 

P 

time  is  then  used  as  the  most  probable  intercollision  time  for  the  Poisson 
distribution  from  which  a  particular  intercollision  time  T(i)  is  chosen. 

Since  the  postcollision  speed  changes  after  each  collision,  the  most  probable 
intercollision  time,  a  function  of  that  speed,  also  changes  after  each 
collision.  Consequently,  after  each  collision,  T(i)  is  chosen  from  a 
different  Poisson  distribution.  In  this  way  the  mild  correlation  between 
intercollision  times  and  postcollision  speeds  found  even  in  the  buffer-gas 
case  is  rigorously  taken  into  account.  Once  the  postcollision  speeds, 
velocities,  and  intercollision  times  are  specified,  x(t)  is  defined  and  may  be 


Fourier-analyzed,  as  was  the  bufferless,  wall-coated  cell,  to  yield  the 
lineshape. 

The  calculation  was  performed  for  a  system  in  which  the  correlation 
should  be  enhanced.  Hydrogen,  with  a  hyperfine  frequency  of  1.42  GHz,  was 
chosen  as  the  emitting  species;  Kr  was  used  as  the  buffer  species.  Because  Kr 
is  larger  in  mass  than  H,  its  buffer  atoms  appear  almost  motionless  to  the 
rapidly  moving  H  atoms.  A  system  temperature  of  25°C  was  used,  and  the 
radiative  lifetime  of  the  H  atom  was  taken  to  be  0.025  s.  The  ratio  of  the 
mean  free  path  to  the  emitted  wavelength  was  set  at  0.01,  placing  the  calcula¬ 
tion  within  the  standard  Dicke  regime.  The  calculated  lineshape  is 
Lorentzian,  with  a  linewidth  of  approximately  220  Hz.  Since  the  linewidth  is 
well  above  that  specified  by  the  atom's  radiative  lifetime,  there  is  no  need 
to  increase  the  radiative  lifetime  employed  in  the  calculation. 

The  buffer-gas  calculation  was  repeated,  but  Instead  of  the  velocity- 
dependent  intercollision  times  based  on  Eqs.  (2.6)  through  (2.8),  speed- 
independent  intercollision  times  were  used.  These  times  were  chosen  randomly 
from  a  Poisson  distribution  whose  most  probable  collision  time  t'  was  given 
by:^9 

(t'  )-1  -  ir  o2N  V  (2.9) 
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and  M*  the  H-Kr  reduced  mass.  This  calculation  also  resulted  in  a  Lorentzian 
lineshape;  however,  the  linewidth  was  nearly  370  Hz.  The  inclusion  of  the 
correlation  effects  yields  a  linewidth  reduced  by  40%  from  the  uncorrelated 
linewidth.  We  find  that  even  in  the  buffer-gas  case  the  effects  of  the 
correlation  between  intercollision  times  and  postcollision  velocities  must  be 
taken  into  account. 


2.3  DISCUSSION  OF  THE  MOTIONAL  NARROWING  PEDESTAL  IN  THE  WALL-COATED  CELL 


As  indicated  in  Section  2.1,  the  trajectory  analysis  we  described  is  in 
excellent  agreement  with  the  observations  of  spectral  line  narrowing,  but 
seems  to  have  some  difficulty  in  exactly  reproducing  the  experimentally 
measured  pedestals  on  which  the  narrowed-  portion  of  the  spectral  line  is 
found.  To  investigate  the  computationally  generated  pedestal,  lineshapes  were 
calculated  for  a  number  of  cells  of  radii  ranging  from  1.5  to  40  cm.  In  Fig. 
5,  several  of  these  curves  are  presented.  As  would  be  expected  from  the  one¬ 
dimensional  analysis,  with  increasing  cell  radius  the  height  of  the  pedestal 
with  respect  to  that  of  the  central  peak  becomes  larger.  In  each  case, 
though,  the  pedestal  remains  non-Gaussian  with  a  width  of  approximately  6  kHz. 

To  ensure  that  this  nonagreement  in  the  pedestal  region  of  the  spectral 
lineshape  was  not  a  computational  artifact,  the  algorithm  was  subjected  to  a 
series  of  tests.  The  linewidth  of  the  central  spike  is  exactly  what  is 
expected  from  the  simple  Fourier  analysis  of  a  sinusoidal  wavetrain  persisting 
for  the  tQ  employed  in  the  calculation.  The  effects  on  the  central  spike  of 
the  frequency  changes  due  to  velocity-changing  collisions  are  completely 
removed.  If  a  single  fixed  atomic  velocity  is  supplied  to  the  Fourier 
transform  algorithm,  spectral  intensity  is  generated  at  the  correct  Doppler- 
shifted  frequency.  As  another  test,  rather  than  generating  velocities  and 
intercollision  times  from  atomic  trajectories  within  the  wall-coated  cell, 
these  quantities  are  chosen  to  be  appropriate  for  a  Rb  atom  surrounded  by  a 
buffer  gas,  as  described  in  Section  2.2.  The  resulting  lineshapes  match 
exactly  the  results  of  the  standard  statistical  analyses  of  Dicke  narrowing. 
The  lineshapes  change  smoothly  from  narrow  Lorentzians  at  smaLl  mean  free 
paths  to  full  Doppler  lineshapes  at  large  mean  free  paths.  In  light  of  these 
tests,  we  feel  that  the  algorithm  is  correctly  calculating  the  raotionally 
narrowed  lineshape  appropriate  to  the  physical  model  we  have  described. 

In  modeling  the  signals  of  interest,  several  simplifying  approximations 
have  also  been  analyzed.  First,  the  radiating  atoms  are  known  to  stick  to  the 
cell  surface  for  times  on  the  order  of  10“^  to  10"  s,  during  which  the 
radiated  frequency  is  significantly  perturbed  by  interactions  with  the  cell 
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walls.  From  the  results  of  Vanier  et  al.^  the  magnitude  of  this  shift  for 
atoms  residing  on  the  surface  for  10~^  s  is  approximately  64  MHz.  To  observe 
the  effects  of  this  process,  calculations  were  performed  in  which,  between 
each  atomic  trajectory,  the  atom  was  assumed  to  reside  on  the  cell  surface  for 
10“*0  s  and  radiate  at  the  shifted  frequency.  Including  this  wall  sticking 
did  not  modify  the  calculated  pedestals. 

An  additional  difference  between  the  calculational  model  and  physical 
reality  is  associated  with  how  the  atom's  loss  of  coherence  is  treated.  The 
coherence  time  specifies  the  width  of  the  central  peak.  In  our  calculation 
the  atom  is  allowed  to  radiate  within  the  cell  for  the  coherence  time  tQ,  and 
all  coherence  is  lost  instantaneously  at  tQ.  In  reality,  coherence  is  lost 
"gradually"  from  small  random  changes  in  the  phase  of  the  radiation  upon  each 
wall  collision.  To  explicitly  include  this  gradual  loss  of  coherence,  an 
additional  random  variable  n(r)  was  included  within  the  argument  of  the  cosine 
in  Eq.  (2.1)  and  was  assumed  to  be  a  Gaussian  random  variable  whose  value 
changed  upon  each  collision  with  the  wall.  With  the  inclusion  of  n(r),  the 
central  peak's  linewidth  could  be  controlled  by  the  magnitude  of  n(f)'s 
variations  rather  than  by  specifying  a  somewhat  artificial  tQ.  However, 
treating  coherence  loss  in  this  manner  does  not  change  the  pedestal. 

Finally,  we  have  Investigated  the  effects  of  nonuniform  microwave 
intensity  across  the  absorption  cell  as  is  the  normal  experimental  condition 
when  horns  are  used  as  part  of  the  microwave  source.  Although  the  analysis 
was  developed  for  a  uniform  strength,  it  is  possible  to  treat  the  variable 
field  strength  in  a  straightforward  manner.  As  an  extreme  case  we  considered 
a  point  source  located  several  centimeters  from  the  storage.  The  resulting 
A-factor  in  Eq.  (2.1)  is  now  position  dependent.  While  Eq.  (2.5)  may  still  be 
broken  into  a  sum  of  integrals,  one  for  each  intercollision  period,  the 
integration  of  each  of  these  integrals  was  carried  out  numerically.  Even  for 
the  extreme  case  of  a  point  microwave  source,  the  resulting  pedestals  were 
Identical  to  those  calculated  previously,  indicating  that  an  incomplete 
description  of  the  microwave  field  strength  variation  within  the  storage  cell 
is  most  likely  not  the  source  of  the  unexpected  pedestal. 
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The  pedestals  obtained  computationally  are  not  well  described  by  either 
Gaussian  or  Lorentzian  lineshapes.  They  are  well  fit,  though,  by  the  sum  of 
two  Gaussians  of  approximately  equal  heights  but  with  differing  widths.  One 
Gaussian  has  a  width  approximately  equal  to  the  expected  Doppler  width  of 
9  kHz;  the  other  has  a  width  of  about  4  kHz.  Only  when  the  cell  radius  is 
significantly  increased  do  the  relative  heights  of  the  two  Gaussians  change. 
For  example,  when  the  cell  radius  is  set  at  40  cm  the  narrow  Gaussian's  height 
is  only  25%  that  of  the  broader  Gaussian. 

An  interesting  observation  concerning  the  pedestal  shape  was  made.  If 
the  atom's  emission  time  and  the  storage  cell  radius  are  properly  chosen,  no 
collision  with  the  wall  will  occur  during  the  emission  period.  Averaging  a 
large  number  of  these  collisionless  trajectories  yields  a  pedestal  of  standard 
Doppler  lineshape.  However,  once  the  trajectories  are  selected  such  that,  on 
the  average,  a  single  wall  collision  occurs  during  each  trajectory,  the  pedes¬ 
tal  reverts  to  the  double  Gaussian  form.  The  fact  that  the  effects  of  colli¬ 
sions  in  the  spherical  cell  become  apparent  after  a  single  collision  is 
consistent  with  the  findings  of  Sykes. ^0  When  studying  motional  narrowing  of 
Zeeman  resonance  lineshapes  in  inhomogeneous  magnetic  fields,  he  found  that 
considering  only  the  initial  motions  of  atoms  within  spherical  cells  up  to 
their  second  wall  collisions  defined  the  lineshapes  quite  adequately. 

At  this  point  we  are  unable  to  explain  the  apparent  discrepancy  between 
the  calculated  pedestal  and  the  experimentally  observed  pedestal.  The  calcu¬ 
lations  accurately  calculate  the  Dicke-narrowed  lineshape.  The  fact  that 
there  are  significant  differences  between  the  calculated  and  observed  pedes¬ 
tals  may  mean  that  the  pedestal  region  of  the  spectral  lineshape  is  extremely 
sensitive  to  the  particular  experimental  conditions.  There  is  a  paucity  of 
detailed  experimental  measurements  of  the  spectral  lineshape 's  pedestal 
region,  since  the  narrowed  portion  of  the  lineshape  has  always  been  of  more 
interest  from  the  aspects  of  both  its  fundamental  physics  and  its  application 
to  various  devices  such  as  atomic  frequency  standards.  Further  experimental 
Investigations  into  the  pedestal  region  of  these  narrowed  spectral  lineshapes 
could  now  prove  extremely  beneficial  in  improving  our  understanding  of  atom- 
atom  and  atom-wall  collisions  and  spectral  line-narrowing  processes. 


3.  SUMMARY 


It  has  been  shown  that  a  single  formalism  can  be  used  to  describe  Dicke 
narrowing,  collisionally  induced  motional  narrowing.  The  extreme  narrowing 
observed  in  the  wall-coated,  bufferless  cell  results  from  the  high  probability 
that  the  precollision  and  postcollision  velocities  will  be  of  opposite  sign, 
thereby  cancelling  Doppler  shifts.  The  pedestals  observed  beneath  the  narrow 
central  spike  are  non-Doppler,  more  accurately  described  by  the  sum  of  two 
Gaussians,  one  with  the  Doppler  width  and  a  second  with  a  width  of  approxi¬ 
mately  half  the  Doppler  width.  The  calculated  form  of  the  pedestal  is  not  in 
agreement  with  the  limited  available  experimental  data  that  show  a  Doppler 
pedestal.  For  atoms  confined  in  a  buffer  gas,  the  effect  of  the  mild  correla¬ 
tion  between  the  postcollision  velocities  and  intercollision  times  has  been 
analyzed.  Including  this  mild  correlation  results  in  a  somewhat  narrower  line 
than  is  predicted  by  the  standard  analyses  of  Dicke  narrowing,  which  neglect 
its  effect. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer’*  for 
national  security  projects,  specializing  In  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  la  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  Individual  laboratories: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy ,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communlcat ions; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


